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Abstract — We analyze a power distribution line with high 
penetration of distributed generation and strong variations of 
power consumption and generation levels. In the presence of 
uncertainty the statistical description of the system is required to 
assess the risks of power outages. In order to find the probability 
of exceeding the constraints for voltage levels we introduce the 
probability distribution of maximal voltage drop and propose 
an algorithm for finding this distribution. The algorithm is 
based on the assumption of random but statistically independent 
distribution of loads on buses. Linear complexity in the number of 
buses is achieved through the dynamic programming technique. 
We illustrate the performance of the algorithm by analyzing a 
simple 4-bus system with high variations of load levels. 

L Introduction 

Ensuring the quality of delivered power is one of the 
main challenges faced by utility companies. Rapid advent of 
renewable generation and electric vehicle technologies will 
inevitably result in a significant increase in the variations 
of power consumption and/or generation Keeping the voltage 
level within the industry constraints will become an even more 
formidable task for the utilities. One of the most challenging 
aspect of this problem is the uncertainty about the load 
structure in the feeder line. Lack of information about the loads 
makes the problem essentially probabilistic and thus requires 
more sophisticated techniques for analysis and control. 

Traditional methods of controling the voltage level on 
distribution feeder lines include but are not limited to line 
regulators and capacitor banks [1]. Whereas the line regulators 
based transformer tap changes are employed for controlling 
the voltage level in the beginning of the line, the switched 
capacitors are usually distributed in the middle of the line, and 
can be used for smoothing out the voltage drop curve via the 
reactive power injections [2], [3]. Penetration of intermittent 
renewable generators and increasing demand in power supply 
will require novel approaches for controlling the voltage level. 
One of the most promising ideas that has been proposed 
recently in this field is the distributed control of reactive power 
flows via the local inverters attached to renewable generators 
[4], [5], [6]. 

In a recent work [7] we have shown that the simple control 
techniques of distributed inverters can be efficient in reducing 
the losses in radial distribution systems. However, the effect 
of inverters on the overall power quality remained an open 
question. One of the main obstacles to the proper assessment 
of the power quality in the system was the lack of theoretical 
methods of voltage drop analysis in the presence of high 



variations of loads. In this paper we attempt first steps of 
solving this problem by introducing a novel algorithm, that 
allows fast prediction of the probability distribution function 
of the maximal voltage drop in the feeder line. Specifically 
the algorithm allows one to transform the distribution of loads 
on individual buses into the distribution of maximal voltage 
drop. 

The structure of this paper is the following: in the section 
[ill introduce the power flow model based on classical 
DistFlow equations. Then, we formally introduce the main 
object this study: the maximal voltage drop. Extension of 
the model to include the stochastic variations of power con- 
sumption/generation are presented in the section [III] We finish 
by presenting the results of simulations and discussing the 
remaining challenges. 

II. Power flow model 

The power flows in the linear distribution line with N buses 
can be described with the DistFlow recurrence equations [2], 
[3]: 

Pk+i=Pk-rk \,2 " -Pk+i (1) 
Qk+i=Qk — Xk ^ qk+1 (2) 
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where Pk, Qk are real and reactive power flows from bus fc — 1 
to k. Vk is the voltage level on bus fc, and pk, qk are the values 
of real and reactive power consumption on bus k. Tk + jxk is 
the complex impendance of the link between the nodes fc and 
fc + 1. These equations have to be solved with two boundary 
conditions: fixed base voltage level at the beginning of the line 
Vq and zero power flux through the virtual link at the end of 
the fine: Pn+i = Qn+i ~ 0. 

For a class of low voltage distribution lines considered 
here the quadratic terms in the equations ([T][3]) are negligibly 
small, so one can use the linearized power flow equations 
where quadratic terms are dropped. Moreover, for the sake 
of simplicity in this work we will restrict the analysis to the 
homogeneous networks where the ratio Xkjvk = a is the same 
for all links. In this case one can combine the values of Pk and 
Qk into a single variable Sk = Pk+ ctQk (and also introduce 
Sfc = Pk + ctQk)- Moreover, as long as the variations of the 
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Fig. 1. Schematic representation of the linear feeder line. Each bus is 
characterized by the voltage level Vfc, consumption of real and reactive 
powers: Ph,Qk- The power flows between the buses are denoted by Pk,Qk- 



voltage are small compared to its base value, one can use 
the approximation Vj^ = Vq^ + 2Vo(Vfe — Vq)- Using these 
approximations one can rewrite the DistFlow equations in the 
following simple form: 



Sk — '5'fc+i + Sfc+i 
Vfe = Vk+i + PkSk- 
where we have introduced — f^/Vo. 



(4) 
(5) 



A. Maximal voltage drop 

The power quality in the line can be characterized by 
the maximal voltage drop level - the difference between the 
base voltage Vq and the minimal voltage in all buses: Aq ~ 
Vq — minfc 14 . In order to find the value of Aq we introduce 
an intermediate object A„ characterizing the maximal voltage 
drop in the feeder line segment n . . . N: 



A„ = Vn - min Vfe 

n<k<N 



(6) 



Note, that A„ is essentially non-negative quantity. The re- 
currence equation for A„ can be easily derived from (j5]) 
by noticing that there are two distinct scenarios. Whenever 
Vk+i+ PkSk > the voltage level on bus k will be larger than 
the minimal one in the segment k . . . N : Vk = Vk+i+PkSk > 
Vk+i - Afe+i = minfe+i<;<jvy/. Therefore, according to 
the definition (|6]l the value of A^ will be given by Afe = 
Afe+i + PkSk- In the second case when Afe+i + pkSk < 
the value of Vk will be smaller than minfe_|_i<;<Ar Vi and thus 
the minimal value of the voltage in the segment k . . . N is 
achieved on bus fc, and thus A^ ~ 0. Formally one can express 
these two scenarios with the following equation: 



Afe = 




if Ak+i + PkSk > 0, 
if A, ■ " - - 



PkSk < 0. 



It is important to note, that although the equations (4]5i for 
Sk and Uk are linear, the resulting equation (|7]i for Afe is 
essentially non-linear because of the second case where Afe 
is reset to zero. This second case is realized only if the total 
power S'fe flowing through the link is negative. This can happen 
only in the presence of buses that inject power in the system. 
As long as 5fe variable incorporates both real and reactive 
power, the injection of S flow can be associated either with 
distributed generators that inject real power or with capacitor 



banks that inject reactive power. In the absence of power 
injection in the system, the voltage is a monotonous function, 
and the problem of finding the maximal voltage drops reduces 
to calculation of Vq — Vn which is a linear function of the 
loads Sfe. This problem can be solved analytically without any 
sophisticated algorithms, proposed in this paper. Here we focus 
on the nontrivial situation, and assume that there are buses that 
inject power in the line. 

By solving the recurrence relations Q and (|7]i backwards 
in k with the intial conditions Sn+i = 0, Ajy+i = one 
can easily find the value of maximal voltage drop Aq in the 
whole network. These equations form a basis for the statistical 
analysis of voltage drop in presence of uncertainties about load 
levels in the system. 

III. Statistical description of power flow 

Whenever the precise values of pk , Qk are not known, it is 
not possible to find the precise value of the maximal voltage 
drop Aq. Instead one has to develop a statistical approach 
to the problem. The traditional way of characterizing the 
uncertainties in the system is to study the probability density 
functions (PDF) of system state. Assuming that the loads 
on different buses change independently one can define the 
probability of observing pk + aqk — s via the PDF 7rfe(s) 

The origin of these distributions is not important for our 
analysis: they could be either prior "Bayesian" distributions 
based on the measurement history, or distributions derived 
from some statistical model of the load on bus k. The statistical 
independce assumption implies that the joint probability of 
observing pi + aqi = s[, . . .pN + aqw — s'^ is given by the 
product 7r(si ... s^) = life Ms'k)- 

The goal of the statistical analysis is to transform the PDF of 
loads into the PDF of the maximal voltage drop Aq i.e. into the 
probability of observing the given value of maximal voltage 
drop. One can easily write the formal integral expression for 
the probability of having Aq < A. This can be done by noting 
that Ao < A if and only if Vfc > Vb - A for all < fc < TV. 
Given the formal solution Vfe = Vfe (si ... sat) of (4 5 1 one can 
formally write: 

Prob(Ao<A)= / lld.Sk7:i.Sk)e{Vk{si....SN)-Vo + A) 
•' k 

(8) 

where d{x) is the Heaviside (Unit step) function. Here and 
throughout the text we assume that intergration is taken 
over the domain (— oo, oo), unless the integration domain is 
specified explicitly. On practical level, this solution is unusable 
for large systems with ^ 1, as long as there are no 
fast "black-box" algorithms to evaluate this multidimensional 
integral in case ^ 1. Even in the simplest case of uniform 
distributions Hkisk) the problem is reduced to calculation of 
the volume of a highly dimensional polytope, for which state 
of the art algorithms require at least 0{N'^) operations [8]. 

'Throughout the paper we use the term probability distiibution to refer to 
the probability density function: by definition the probability of observing 
Sfe £ [s — "5/2, s + <5/2] is equal to 7rj;(s)5 in the limit of 5 — > 



It is therefore important to develop alternative approaches 
to computing the PDF of Aq that would require only linear 
in N number of operations and would be applicable even for 
large systems. We propose a specific algorithm of the kind 
that exploits the radial structure of the network and statistical 
independence of load levels on different buses. The main idea 
behind the algorithm is to transform the recurrence relations 
( 4]7 1 into the corresponding relations for the joint probability 
functions of Ak,Sk and solve them in an iterative way. 
This general strategy of reusing the previous computations 
via decomposing the problem in recurrence type relations is 
usually referred as dynamic programming [9]. 

Existence of recurrence relations (4]7i for the pair of 
variables Sk,Ak which expresses the PDF 11^(5', A) of the 
values of S'a , A^ via the PDF 11^+1(5', A) of the values of 
Sk+i, Afe+i. This Chapman-Kolmogorov type relation can be 
written in a most compact way with the use of Dirac Delta 
Functions S{x): 
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Fig. 2. Probability distribution function of the maximal voltage drop Aq. 
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[e{A' + pkS)d{A - A' - pkS) + e{-A' - pkS)5{A)] (9) ^^N+ii^) = 0. In the next section we will describe the results 



where we have used 5', A' to denote the values of Sk+i, Ak+i 
respectively. This recurrence relation has to be solved with 
the initial condition nAr+i(S', A) = 5{S)5{A). Formally 
this equation completes the construction of the algorithm, as 
one can iterate it N times to obtain the PDF of maximal 
voltage drop in the system Aq. For a given target precision, 
the complexity of performing single iteration depends only 
on the numerical discretization of 11^(5', A) and does not 
depend on the total number of buses in the system. Therefore 
the total number of iterations grows like 0{N). However, 
implementation of a single iteration (j9]l can be obstructed by 
a non-analytic nature of the resulting joint PDF life (5, A). 
Performing the first few iterations analytically, one can show 
that for continous distributions 7rfe(s) the general expression 
for life (5*, A) can be decomposed in three parts: 

Hfe {S, A) = {S, A) + n W (5),5( A) + n(') {S)5{A - 2rkS) 

(10) 

where J\%{S,A) is analytic part of the PDF, and are 
the prefactors in front of the non-analytic ones. Substituting 
this decomposition in (j9|l one obtains the coupled system of 
equations for n%{S,A) and n[^^'^^ 

ni{S,A)^e{A) J dS'7rk+i{S^S')Ul^,{S',A-2rkS) + 
e(A)9{A-S) A-2rfe5 



of solving this system for a sample 4 bus system. 

IV. Results 

Efficient implementation of the proposed algorithm suit- 
able for large heterogeneous systems is a challenging task 
that requires a thoughtful selection of suitable discretization 
technique for the PDFs, that would allow fast evaluation of 
convolution integrals in ( [TT] i and adaptation of the discretiza- 
tion domain to the width of the distribution functions that will 
grow with each iteration. Analysis of possible approaches to 
these problems is beyond the scope of this paper. Instead, here 
we report our results of a "proof of concept" study of toy 
feeder line system consisting of iV = 4 buses where for the 
sake of simplicity we have set all pk to pfc = ViT^-p.u.jkW 
and the load distribution functions to 



[exp(l.Os) if s < 0. 
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Note, that this distribution function implies that flux of real 
and reactive power can have an arbitrary sign. Although the 
probability of power consumption (s > 0) is higher, there is 
also a finite probability of a given bus injecting the power 
in the line. The average value of power consumption was 
set to s = 2.0kW whereas the standard deviation of power 
consumption is equal to (ig w S.WkW, so the fluctuations of 
power consumption/generation are indeed strong in the system. 

In order to avoid technical difficulties associated with dis- 
cretization of the probability distribution functions, we have 
solved the equations (111 analytically using the Wolfram 



n 



(S) = e{s) I ds'TTk+i{s-s')ni%{s') (ii) 



Mathematica computer algebra system. On each step the 
functions 11^.(5, A), 11^. ' {S) were represented as a piecewise 
analytical functions, and the convolutions ( [TT] l were performed 
rigorously without any numerical approximations. Although 
the number of terms in the analytical expressions was growing 
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Fig. 3. Contour plot of the joint probability distribution nQ(5, A) of the 
total load So and maximal voltage drop Aq. 

quite rapidly, it was possible to find the expressions for 
ng(S',A) and n[,^'^^(5') for a system with N = 4 buses. 
These expressions were used to calculate several statistical 
characteristics of the system. 

The probability of exceeding the voltage level constraints 
can be calculated from the analytical part of the probability 
distribution of the real voltage drop Ao or equivalently the 
probability density in the region of positive voltage drops: 

P(Ao|Ao >0) -n{,2^(Ao) + I dS'U-oiS',^o) (13) 

This corresponding result is presented on the figure |2] As one 
can see the distribution is centered around characteristic value 
of Aq w A.lp.u., however the variations of the typical voltage 
drop are very high, for instance the probability of exceeding 
twice the average value Aq > 2Aq is about 13%. This result 
confirms the central thesis of the study: variations in load 
levels result in large variations of the voltage drop along the 
line and thus significantly increase the probability of exceeding 
the limits set by regulators. 

Another interesting object is the joint distribution of the total 
power consumption in the system 5*0 and the maximal voltage 
drop Aq. The corresponding distribution nQ(S'o, Aq) is shown 
on the figure [3] As one can see there is a strong correlation 
between 5*0 and Aq. This correlation can be potentially used 
for designing the voltage control techniques that respond to 
the total load So and adjust the voltage Vq to suppress the 
risk of exceeding the allowev voltage drop limits. Statistical 
analysis of various control techniques is an interesting problem 
that can be approached with the techniques introduced in this 
article. 



V. Discussion 

Algorithm presented in this paper is a first step on a 
long path of development a new generation of tools that 
would allow the power system designers and controllers to 
assess the risks associated with fluctuations of load levels 
and uncertainties in the system. Fast ways of calculating 
the probabilities of high voltage drops can be utilized in 
system design optimization packages or stochastic control 
processes. One of the specific applications of this algorithm 
could be the analysis of the effect of distributed renewable 
generation on the power quality in the distribution system. The 
probability distribution functions of voltage drop can be used 
for determination of the critical penetration levels of renewable 
generators that can be sustained by the current lines. It can be 
also used for comparison of different techniques of voltage 
control with reactive power generators: capacitor banks [2] or 
small-scale inverters attached to photovoltaic elements [7]. 

However, there are still several issues that have to be 
dealt with before the algorithm can be put to use. The most 
important of them is the choice of suitable discretization 
technique, that would allow efficient evaluation of integrals 
( [TT| ). The main obstacle here is the piecewise-analytic struc- 
ture of the convolution kernels that will result in piecewise 
continous structure of the distribution functions. This structure 
suggests that the probability distributions should be discretized 
with finite element partitioning of the domain. However, the 
composition of domains will have to be adaptable, as the 
structure of the distribution function will change with each 
iteration. For practical problems it is also important to extend 
the current algorithm to analyze the joint statistics not only 
of the minimal but also the maximal levels of voltage along 
the line. Although on formal mathematical level this extension 
is straightforward, the numerical realization will become more 
complicated, as the distribution functions will depend on three 
variables: power flow S and maximal and minimal voltage 
drops. 

VI. Conclusion 

We have presented a dynamic programming algorithm for 
calculating the probability distribution function of the maximal 
voltage drop in the linear feeder line. The complexity of 
the algorithm scales linear with the number of buses in the 
system, which makes it superior in comparison to the standard 
black-box type algorithm for solving similar problems. The 
proposed algorithm requires iterative application of integral 
operators that relate the probability distributions of local 
voltage drops on neighboring nodes. Efficient implementation 
of this convolution is a difficult problem that will be hopefully 
solved in forthcoming works. In this paper we have tested the 
algorithm on a toy 4-bus homogeneous system. The resulting 
distribution functions show that the variations in load lead to 
variations in voltage drop. Moreover they also show strong 
correlation between the voltage drop and total load of the 
system. 
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